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Abstract 

In this paper, we examine interwell tunneling between a pair of fractional 
quantum Hall liquids in a double quantum well system in a tilted magnetic 
field. Using a variational Monte Carlo method, we calculate moments of the 
intra-Landau level tunneling spectrum as a function of in-plane field compo- 
nent Bu and interwell spacing d. This is done for variety of incompressible 
states including a pair of v = 1/3 layers ([330]), pair of v = 1/5 layers ([550]), 
and Halperin's [331] state. The results suggest a technique to extract interwell 
correlations from the tunneling spectral data. 

PACS: 73.40Gk, 73.40.Hm, 73.20.Dx 
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I. INTRODUCTION 



Traditionally experimental studies of the quantum Hall effect have been restricted to 
magnetotransport, surface acoustic wave, and capacitance spectroscopy studies. However 
during the last several years new spectroscopic techniques have been developed. These tech- 
niques include photoluminescence, inelastic light scattering!!!, and tunneling spectroscopy^. 
The latter has just begun to be used to investigate the tunneling spectra between a pair of 
weakly coupled two-dimensional electron liquids!'!!! and between a two-dimensional liquid 
and a three dimensional doped substrateS. 

This paper will be primarily concerned with tunneling between a pair of parallel quantum 
wells in the high field regime. Experiments of this sort have been performed by Eisenstein, 
Pfeiffer, and West!. For filling fractions ranging from v = 0.48 to v — 0.83 per layer, 
these authors find the following features in the I-V characteristics. At low voltages, the I-V 
characteristics exhibit a pseudogap with activated temperature dependence. (The activation 
temperature is about 5-10K.) Above this pseudo-gap, Eisenstein et all find a featureless band 
of intra-Landau level excitations which peaks at a voltage 2A2 ~ 0.45e 2 /e/ m . At still higher 
voltages, they find a second gap which separates the intra-Landau excitations from the inter- 
Landau level and inter-subband excitations. In view of the above experimental results, it 
is clear that a detailed theoretical understanding of inter-well tunneling in the high field 
regime would be desirable. 

To date, the theory of inter-well tunneling in the high field regime has focussed on 
dynamical issues like the size and origin of the pseudo-gapii. For instance, Johannson and 
Kinaret! found a tunneling pseudogap in the non-linear I-V characteristic using a Wigner 
crystal model. Above the pseudo-gap, they find a featureless band similar to that found 
in the experiment. In addition to this work, there are a number of exact diagonalization 
calculations of the single electron spectral functions!'! which may be used to calculate the 
tunneling conductance. Because of finite size effects, the diagonalization calculations do 
not obtain continuous I-V curves. Nevertheless, these calculations and the Wigner crystal 
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model all seem to obtain a peak in the tunnel current at a voltage which is consistent with 
experiment. In addition Yang and MacDonaldS have calculated the tunneling density of 
states of a disordered 2D electron gas in a strong magnetic field. They find a suppression 
of the tunneling conductance at small voltages but no pseudogap. He, Platzman, and 
Halperini also find pseudo-gap behavior in a pair of v = 1/2 Halperin Lee Read0 Fermi 
liquids. Finally, Efros and PikusS have studied a lattice-gas model of a classical electron 
liquid using a Monte-Carlo methods. 

In contrast to those references&~0 which focus on dynamical issues, this paper will fo- 
cus on features of the equilibrium interwell correlations which may be extracted from the 
tunneling conductance. We will argue that a significant understanding of static correlations 
may be achieved by studying the tilt angle dependence of a few spectral moments extracted 
from experiments using a tilted field geometry. The case for this will be made through an 
examination of tunneling spectral moments obtained from a variational Monte Carlo (vMC) 
calculation. The vMC method described here, we feel, gives results which are complementary 
to those obtained from exact diagonalization and independent boson model calculations. Of 
course the vMC method does suffer from certain limitations. Most notably vMC requires 
the use of reasonably manageable wavef unctions. This requirement will restrict our investi- 
gations to the simplest incompressible states which can occur in double well systems. One 
might be concerned that a comparison of these results with the experimental data would be 
meaningless since the data is restricted to compressible states whose filling fraction ranges 
from .48 to .83 per layer. We would argue however that this is not a serious problem for two 
reasons: First, the experimental data of Eisenstein et all is found to be rather insensitive 
to the filling fraction. Secondly, the results reported below involve the first few spectral 
moments which we will argue are primarily sensitive to short range interwell correlations 
and are less sensitive to to slight changes in filling fraction. 

The need to work with simple trial wavefunctions gives rise to a second problem which 
we will now describe: Consider the tunneling between a pair of v = 1/3 states. A question 
one might like to ask is what sort of inter-well correlations will be induced and how might 



one detect them in a tunneling experiment. The difficulty which is encountered is that 
the wavefunction which includes the correlations is not simply the product of two v — 1/3 
Laughlin wavefunction, it is some perturbation of the product wavefunction. Unfortunately, 
identifying the form of the perturbed wavefunction and using it in a calculation would add 
significantly to the complexity of the discussion. To avoid this problem we proceed as follows: 
(i.) First, we will calculate the relevant spectral moments using unperturbed wavefunctions 
like a pair of v = 1/3 [330] and a pair of v = 1/5 [550]. The results of this calculation 
may then be compared with experiment in order to address issues regarding the quality 
of data and the quality of the theoretical model, (ii.) Secondly we address the issue of 
the identification and characterization of inter-well correlations by studying the HalperinEl 
[331] state. This wavefunction displays significant inter-well correlations and the required 
calculations are straightforward to implement. The results of the calculation on the [331] 
state suggests a general approach to analysing tunneling data. We believe the approach 
should be able to identify inter-well correlations independent of whether they are present in 
the zeroth order wavefunction (as is the case with the [331]) or are present only in the exact 
wavefunction (as is the case with a pair of v = 1/3 liquids). Of course the calculation of 
spectral moments obtained for the [331] state could also be compared to experimental data, 
should such data ever become available. 

Consider, therefore, the device shown in fig. 1 on which a tilted magnetic field B = 
(B\\,0, B^) is applied. We will model the double well system using the Hamiltonian H = 
H Q + H t where H is defined to be 

r _, g \ r 

H = Yl / *"c4(V - i— A) 2 c a + - V / dr 1 dr 2 V a p(r 1 - r 2 )4(ri)cl(r 2 )c / 3(r 2 )c a (ri) 

Zm eff a J hc 1 a fi J 

(1) 



where H t = J d 2 r[t exp(z^ / dz A z (r)) ctcj + h.c] is the tunneling term. The electron- 
electron interaction is taken to be V^(r) = Vn(r) = e 2 /er and V^j_(r) = e 2 /e(r 2 + d 2 ) 1 ^ 2 
where d is the interwell spacing and e is the dielectric constant. 

Our discussions will focus on pairs of v = 1/m states (also denoted [mmO]), and the 



Halperin [331] state recently observed by Eisenstein et allI2l in double well systems and by 
Suen et alu in wide single quantum wells. The Halperin [mmn 1 statesMzl are described by 
the wavef unction 

V mrnn = H( Zi - Zj ) m ^ I] e-!^ |2/4 ^ (2) 

i<j i 

where 

{m <Tj = Oj 
(3) 
n a-i^ o, } 

and where a is a pseudospin index which labels the two wells. For simplicity, we will assume 
that the actual electron spin is frozen out by the magnetic field. The [mmn] wavefunction 
describes an incompressible state with a total filling factor v = i/f = 27rl 2 n nQ = 2/(m + n) 
where n is the total electron density on the two layers and where l m = (hc/eB±) 1 / 2 . 

As a function of the voltage bias V, between the two wells, we wish to calculate the 
inter- well tunneling current I t . To do this we first define the tunneling operator 

S+(q)= J d 2 r c+ q exp dB\\ y^j (4) 

where 

^ = 2udB\\/%y 

Expanded in powers of to, the tunneling current is 

It = e< j z >v + 2e |t | 2 Im[X ret (eV )} + ... (5) 

where 

/oo 
dte^i9(t) ([S+(q,t), S+(q,0)}) (6) 
-oo 

In eqn. ||, the first term would describe a Josephson effect, if such exists, the second term 
describes incoherent tunneling. In the following discussion, we will only consider unpolarized 
states m ^ n where the first term in eqn. |^ vanishes. In this case, no Josephson effect effect 
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will occur. (We refer the reader to appendix B for a more detailed discussion regarding the 
absence of a Josephson effect.) 

Now using variational Monte Carlo, we will calculate various spectral moments of the 
form 



W k (B ]{ ) = J d(eV) I t {V,B{){eV) k 



(7) 



The k = moment is W (-B||) = 27re|io| 2 -/VS'++(27r£?||<i/0o) where the pseudospin correlation 
function is defined by 



$±± (/>) = jj (sl(k) S±(k) 



(8) 



where S+(q) is the Fourier transform of the interwell tunneling operator, etc;, projected 
onto the lowest Landau level. The k — 1 moment is Wi(B\\) = 2ire\to\ 2 N f ++ (2irB\\d/(j)o) 
where the oscillator strength is 



2N 



Si (q),\H, S + (q) 



(9) 



and where N is the total number of electrons in either well. Results for Wi(B\\), will 
be presented in terms of < eV > the mean voltage bias which is defined by < eV >= 
Wi(B\\) /Wo(B\\) = f ++ (q)/s ++ (q). This expression is exact to the extent that the \mmn > 
wavefunction is the exact ground state. < eV > may be written in the more explicit form 



< eV >-- 



Si (q), \H, S+(q) 



Sl(q) S±(q) 



(10) 



This expression is, of course, the basis of the single mode approximation]!! We do not, 
however, refer to the variational Monte Carlo calculation as a single mode approximation 
since it does not assume a single collective mode. 

The oscillator strength f++(q) may be calculated using the following expression: 

1 r d 2 k 



/++(*) 



{a a p (k, q) S a/3 (k) + b(k-q,q) [S ++ (k) + (k)] } (11) 



4 J (2nY 

where the various quantities will be defined as follows: S a p(q) is a structure factor matrix 
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S a /3(k) = — (n a (-k)np (k) 



and b(q, k) and a a/ g(fc, q) are 

6(ife,g) = 2V n (k)e {kq + kq)/2 - V n (k) (e q k + e k q ) 

and 

a af 3 {k, q) = Var, ( k ) m iP ( fc > <?) + m "7 ?) 

where 



(12) 



(13) 



(14) 



m(fc, g) 



1 — exp (kq — qk)/2 
c.c. 1 



exp(-|g| 2 /2). 



(15) 



This expression is obtained from the definition [equation (§)] using manipulations of the sort 



described in reference 18. 



The structure factor matrix is been calculated for the [330] and [550] states in ref. [18|. It 



has also been calculated for the [331] state in ref. [19|. So only the calculation of S ++ (k) = 

S (k) will need to be done. This is the task of the next section. Readers uninterested in 

the technical details of this may simply skip that section. 



II. MONTE CARLO CALCULATION OF THE PSEUDO-SPIN CORRELATOR 



This section describes the calculation of S ++ (k) using the variational Monte Carlo 
method. To do this we first introduce the pair spin correlation function: 



g++{r) = tt 



< Ni,N 2 


S 1 _(z')S 2 + (z)\N 1 N 2 > 


< N^z 


NxN 2 > 



(16) 



where S ± (r) denotes the spin of the i-th electron projected to the lowest Landau level, i.e. 



S± = a l ±5(r — Ti) where 5{r — r«) is the Fourier transform of e Sz j e z k z i . The pair spin 



correlation function is related to S ++ (k) as follows: 



S++ (k) = x^' 



\k\ 2 12 



+ noX^Xi Id r e 



-ik ■ r 



9++ [r 



(17) 



In equation (|T6| ) Q denotes the total area of the two dimensional electron gas, and x a is the 
fraction of electrons on the a well. 

Using the explicit form of ^ m mn and assuming balanced quantum wells, i.e. — iVj_, 
one can write the pair spin correlation function in the form 



9++{ z ~ z 



\z' - z^ 2 



j^z^ < 5^ ai 5^ 2 5(z ! - z x )8{z - z 2 )A zlZ2 [ Zi ] > N 



(18) 



where the weight factor A Z1Z2 [z k ] is defined by 



A 



ZlZ2 



n (Zi - Z k ) (Z 2 - Z k ) 
k=N ] +2 

n {zi - z k ) (z 2 - z k ) 

k=3 



(19) 



and where < (...) >n denotes 

(iV T + l,iV i -l|(...)|iV T + 1,^-1) 



(20) 



(JV T + 1, N t - 1|JV T + 1, JVj, - 1) 

The origin of the factor A ZlZ2 [z k ] in equation ( |18D may be understood as follows: First 
S + (z) \Ni,N±) 7^ |JV| + 1, N± — 1) because the electron transfer from bottom to top well 
leaves m — n more zeros than the n required for a \mmn > state. Moreover, the new 
electron on the top well is bound to n zeros rather than the m zeros characteristic of a 
\mmn > state. As a result one must multiply, \N^ + 1, iV| — 1 > by an additional Jastrow 
factor in order to obtain S+(z)\Nj, N± >. This extra Jastrow factor leads immediately to 
factor A zlZ2 [z k ]. The manipulations which lead to the complete expression presented in 
equation (fl8|) may be found in appendix A. 

The expectation on the right hand side of equation (|18D may be obtained from a simu- 
lation of a pair of mobile impurities in a two component background plasma. To do this we 
rewrite the equation as 

A zlZ2 [zi] \ 

I ~ ~ 1 2(m— n) I 
\Z\ - Z 2 \ v IN 



9i 



\z - z') = n 2 (-l) n ~ m U ai 5 1a2 5(z' - z x )8{z - Z2 )- 



(21) 



Next, we split A ZlZ2 [z^ into its modulus and a part with modulus 1: 
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A ZlZ2 [Zk] — A ZlZ2 [zk] 



Nf+Ni 

II \zi — Zf~\ \Z% — Zk\ 
k=N r +2 

II \Z\ — zA \z 2 — Zk\ 



(22) 



fc=3 

Then, using importance sampling, we absorb the second factor in equation (|2"2"| ) as well as 
the factor of l/|zi — z 2 \ 2( - m ~ n ' > into a new analogue plasma Hamiltonian for the Monte Carlo, 
leaving only A ZlZ2 [ z k] to average over. The classical Hamiltonian of the (modified) analogue 
plasma is 

PH plasma = E Qij In Tij + E r?/2 (23) 

i<j i 

The interaction strength between a pair of background ions is 2m(o"j, <jj). For an impurity 
and a background ion Qu = Q 2 i = m + n. Between the two impurities Q12 — 2n. 
During the simulation we calculate the pair spin correlation function using 

9++(r n ) = 1 — — ■ B ++ (r n ) (24) 

A z ir[n z — (n — ±) z \x z noN* 

where the new notation is defined as follows: 

N MC 2 ' 

N* - E E 1 (25) 
fc=i i=i 

where i = 1, 2 labels the impurity ions and k labels the sampled configurations. The prime 
indicates that only configurations where the impurity ions lie in a circle of radius i?* centered 
at the origin are to be included. B ++ (r n ) is a bin counter which keep track of the contribution 
to the sum of A ZlZ2 [zi(k)} associated with impurity pairs with separation z% — z 2 such that 
nA < \zi — z 2 \ < (n + 1)A where A is the bin width. More precisely, B ++ (r n ) is defined by 

Nmc ' „ 

B ++ (r n ) = E A ZlZ2 [z t (k)]6(\ Zl -z 2 \- nA)6((n + 1)A - \ Zl - z 2 \). (26) 
fc=i 

where Zi(k) denotes the position of the i-th ion in the k-th sampled Monte Carlo cycle. 

We used the Metropolis algorithm to calculate B ++ {r n ). During each cycle, we attempt 
as many Monte Carlo moves as there are particles. We began our simulation with an initial 
equilibration period of 10 3 cycles. During this equilibration period, the ion step size was 



adjusted until an average acceptance ratio of 0.5 was reached. The step size was kept fixed 
after the equilibration period. We then ran for another 2 x 10 6 cycles, sampling one out of 
every ten cycles. This gave a total of Nmc = 2 x 10 5 . Usually in simulations of this sort, 
one counts only those pairs where one of the impurities lies in a circle of radius R*. In this 
way one can reduce finite size corrections to the pair spin correlation function.0 Various 
choices for R*, system sizes, run lengths were tried. Ultimately, we concluded that R* — > oo 
on a system of 200 particles per well gave an acceptable balance of systematic finite size 
errors and statistical error. Specifically, we concluded that the finite size error is smaller 
than the noise due to the Monte Carlo procedure. The statistical error in g ++ (r) for the run 
of 2 x 10 6 cycles was 0.003 for < r/l m < 4.5 except near r/l m = 1.7 where it increases to 
0.008. These represent percentage errors of 1-3%. 

The results for [330] and [331] are presented in figure 0. First consider [330] data. In the 
absence of interlayer correlations one can calculate g ++ (r) exactly. To do this, one notes!! 
that for any pair of uncorrelated liquid states 

< 84z)S + {z>) > = < [ct(z)c + (z)][cl(z')c4z')} > 

= -^)(2Tr)~ 2 exp-\z- z'\ 2 /2 (27) 



This is obtained by Wick factorizing the left hand side into products of the single particle 
matrix. Such a manipulation is exact in the absence of interwell correlations. One then uses 
this result together with eqn. |17] to find that 

g ++ (r) = -e- 2 / 2 (28) 

This result is the solid line through the [330] data presented in fig. 0. The agreement 
between Monte Carlo data and the analytical theory is quite satisfactory. 

Next consider the [331] Monte Carlo data. As is the case with [330], the pair spin 
correlation function for [331] is short ranged. It peaks at r = 1.7/ m and is quite small 
for r > 5.0/ m . This behavior is quite different than what would be obtained from an 
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[mmm] state. In particular if n — > m, then A zz /[zk] — > 1 and lim g ++ (r) — > x^. This 
demonstrates that [mmm] exhibits ODLRO. The presence of ODLRO in the [mmm] state 
has been discussed by Zee and WenH in the contex of a possible Josephson effect for this 
system. Given the absence of the ODLRO in the [331], we conclude that a Wen-Zee type 
Josephson effect will not occur for this state. 

Instead of working directly with the numerical g ++ (r), it is convenient to work with an 
analytic fit. A convenient choice to fit g ++ (r) is 

1 



g++{r) = a(r + br ) exp 



(29) 



To within the accuracy of the Monte Carlo data, one can fit g ++ (r) from the simulation 
with 200 particles per well with a = 0.130, b = 0.0 and s = 1.189. The fit for g++(r) is the 
solid line passing through the [331] data in fig. ||. 

Having obtained the analytic fit for g ++ (r), we insert this into equation (|T7| ) to obtain 
the spin correlation function S ++ (k) and the integrated spectral weight W (B^). See fig. |3|. 
Of course, for the [mmO] state Wo may be obtained analytically using eqn. |2l| In fact, one 
can use eqn. [28| to obtain Wq(B\\) for any pair of parallel liquid states with filling fractions 
{y\,vi) provided that interwell correlations can be ignored. The un-correlated limit of Wo 
will be denoted by W nc where 

W nc (B\\) e\t \ 2 , , 1 d B\\ 2 

~^r~ = nii^ 1 ~ ^ )exp " [ 2 1 5i } ] (30) 

W nc is also plotted in fig. [| The results for W (5||) will be discussed in detail in the next 
section. 



III. IN-PLANE FIELD DEPENDENCE OF THE TUNNELING SPECTRUM 

In figure |3|, we present the integrated spectral weight Wo(B\\) which was obtained from 
the Monte Carlo procedure described in the previous section. We observe that interwell 
tunneling is more suppressed for the [330] than for the [331]. This is due to the strong 
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interlayer correlations in the [331] state: The correlation hole which occurs in the opposite 
layer increases the number of unblocked final states. To understand this, consider the 
following argument: Momentum conservation requires that electrons tunnel along the tilted 
magnetic field lines. Hence even in the absence of manybody correlations, tunneling will be 
suppressed because the relevant matrix element involves the overlap of a pair of Gaussian 
wavefunctions displaced by d* = d(B±/B\\). See fig. £|. Thus eqn. |3C] is telling us that 
for a pair v = l/m states, the integrated tunneling conductance is determine by this single 
electron effect. Intrawell correlations are irrelevant to Wq since the only approximation 
which went into the derivation of W nc was the Wick factorization of the spin correlator in 
eqn. 



30. 



Contrast this state of affairs with what happens in the [331] state (or any other [mmn] 
state with n^O). In this case, tunneling will be enhanced when B\\ = 0, because the electron 
will tunnel directly into its correlation hole on the opposite quantum well. However, if the 
field is tilted, the electron will miss the correlation hole. Hence an in-plane magnetic field 
will reduce Wq. See fig. [5[ To separate correlation effects from the suppression of single 
electron tunneling, we define the 'correlation enhancement ratio'i? = Wq(B\\) /W nc (B\\) . In 
fig. [], we plot R vs B\\. Also plotted in fig. ^| is the radial distribution function gn(d*) 
which is defined by 



x a X/3n g a p(r) = — / ^ ~ ( 31 ) 



where a, (3 =|, j. This result for g-\i(d*) was previously obtained in ref. [L^. We see from 
fig. [5] that the correlation enhancement is maximum when d* lies within the correlation hole 
and has a minimum when d* is near 3.2l m , the radius of the first coordination shell. 

So far the discussion of identifying the correlation hole by using Wq(B\\) has focused 
on the [331] state. However, the discussion is in fact a bit more general: Indeed, even 
though the [mmO] wavefunction gives R — 1, the exact ground state would deviate from the 
[mmO] state. In this case@, will develop a weak correlation hole which deepens as d/l m 
decreases. Presumably R would then reflect the existence of the induced correlation hole. 
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Next we wish to consider < eV >. To calculate < eV > we first evaluate f++(q) using 



eqn. 11 and then we use eqn. [TO to obtain < eV >. Numerical quadratures were verified by 



comparing them with analytical results valid in the d = limit. The sensitivity of < eV > 



to the choice of fit parameters (a, &, s) (see eqn. ^9[) was also studied. Typically the error in 
< eV > due to different choices of fit parameters was less than 5% for changes of up to 10% 
in the fitting parameters. The best fit results are presented in figs. ^ and |[a-c for various 
well spacings for the [330], [550] and [331] states. 

The first thing that one notices in fig. [8] is the strong dependence of < eV > on the 
layer spacing: As d — > 0, < eV > collapses when B\\ = 0. The reason for this is simple. 
For d = and to = the Hamiltonian is SU(2) invariant but the [331] wavefunction is not. 
As a result 5 + (g)|331 > would be a Goldstone mode which implies that < eV > vanishes 
at q = 0. Note however, that numerical diagonalizationJiH indicate that the [331] state 
becomes unstable for d/l m less than some critical value around 0.5. The likely explanation 
is that the band edge of /(V), i.e. 2Ai collapses before < eV > does. 

Next we observe that the B» dependence of < eV > is more rapid for the [331] state than 
for the [330] state. This rapid dependence of < eV > on B» coincides with d* moving out of 
the correlation hole. Once cf* moves beyond the first coordination shell interwell correlations 
become irrelevant and the dependence of < eV > on Bu will be similar to that obtained for 
a pair of uncorrelated fluids i.e. like [330]. This then explains the similarity of the < eV > 
vs. Bn curves obtained for [331] and [330] (see fig. 0) when ql m > 3. 

To summarize, we have computed tunneling spectral moments W and W\ of the [330], 
[550], and [331] states as a function of the in-plane magnetic field in a tilted field geometry. 
We argue that the ratio R = Wo(B\\) /W nc (B\\) provides a qualitative method for imaging 
the correlation hole in g^(d*). We have also studied the behavior of < eV >= Wi/W . The 
results presented in fig. |8]show that the mean of the intra-Landau level spectrum will rapidly 
increase as = dB\\/B± moves through the first configuration shell. Finally, arguments 
presented in appendix B demonstrate that no Josephson effect of the sort proposed by Wen 
and Zee can occur if m ^ n. 
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APPENDIX: A. DERIVATION OF PAIR SPIN CORRELATION FUNCTION 



In this appendix, we will derive the explicit form of the pair spin correlation function, 



i.e. equation (18) from its definition given in eqn. [16|. To recast eqn. [16| into the desired 



form, we will need the identity 

S+iz)^,^ >= (-l)"A<«+D [nj(* - Zj ) A ^p\(z)} \N t + l,N l -l> (Al) 

where Aj = m(|, <jj) — m(f, cr,) and where U'j means to omit the j = i factor. The systematic 
derivation which gives the above result is straightforward but will not be given, since a 
few minutes reflection should convince the reader that the above expression is correct. In 
particular, the Jastrow factor on the left side of equation ( |A1| ) has already been discussed, 
so one only needs to consider the factor p\(z). This factor ensures that the state on the 
right hand side has the z th particle located in the top well at position z which, of course, 
is the case for S + (z)\ N^, iVj_ >. 

Using this identity one then obtains 

< N^NzlS^slWlNxNz >= 

< AT T + 1, Nl - l|ni(2* - z^f^liz)^ - z t )^ |JV T + 1, N { - 1 > (A2) 

In the above expression, one may replace with for n = 1,2. Then, after some straigh- 
forward manipulations of equation (|X2j), one may write g++(r) in the form 
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Q-*g ++ (z> - z) = ^_ l) "?) [I 3 N] < 5 l 4 2 S - zi) 5 (z - z 2 ) AM > N (A3) 



where Q[N a ] =< A^.|A^ > is the normalization of the Halperin wavefunction, (iV£) = 
(iV| + 1, JV; — 1), and where A zz /[Zk\ was defined in equation (|l9|) . 

According to the plasma analogy, one may interpret Q[N a ] as the configuration integral of 
the two component impurity-free plasma. Q[N a ] is related to the classical partition function 
by 

where A is the (arbitrarily chosen) thermal wavelength of the analog plasma. Using this 
expression one can write the ratio of configuration integrals which appears in equation (|A3| ) 

as 

Q[N 1 + 1, Nl -1] ^ <n) 

This is valid assuming that N a » 1. 

Next we insert the configuration integral ratio into equation ( |A3|) to get the general 
expression for the pair spin correlation function: 

/ -, \n—m t>t 

n- 2 g ++ (z 1 - z) = l ~ -I e+K^) < 5 ^ _ ^ 5 (z _ ^ A >N 

\z' — z\ y ' Ml 

(A6) 

This is valid even if the wells are out of balance. For problems in which the wells are in 
balance, /if = /ij and the above expression simplifies to equation ([TJj). 



APPENDIX: B. ABSENCE OF THE JOSEPHSON EFFECT WHEN M^N 

Recently Wen and Zee have suggested that, if one could separately contact the two wells, 
then at T = the [mmm] states might exhibit a Josephson effect. Unfortunately, it seems 
unlikely that the Josephson effect exists at finite temperature since thermal fluctuations 
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would cause < S±(z) > to vanish and < S + (z)S-(z') > to decay algebraically. More likely 
would be some sort of fluctuation contribution to the tunneling current!!. 

Nevertheless, it is interesting to consider whether a similar T=0 Josephson effect should 
exist in a more general [mmn] state. Below we present a simple argument which strongly 
suggests that the T=0 Josephson effect does not exist unless m = n. 

The argument begins by constructing a simple tunnelling Hamiltonian valid for the sub- 
space of \mmn > states of the form \N^,Ni >. This ignores high energy bulk excitations 
and all edge excitations with non-vanishing wavevector. Next we define a new basis 

\ e h°l>= E a ( Ar i> Ar 2) expi(JV T 0i + N l 6 2 )\N h N l > (Bl) 

Ni,N 2 

where a(N^, N±) is peaked about < Aj > and < ATj_ > whose width obeys N a >> AN a >> 1. 
The tunneling Hamiltonian is 

H eff = It j cos(fl T - ; ) - /i T iV T - ^ [ N i (B2) 

where tj =< iVj + 1, — l\H t \N^, N± > and where 8 a and N a obey canonical commutation 
relations. . 

According to this rather simple formalism one expects the various Josephson effects if 
tj 7^ 0. To see if this is the case we compare the angular momentum of |A^,A^ > and 

+ 1, Ni — 1 >. One readily finds that angular momentum difference between these two 
states is 

AL Z = (m-n)(Ni - + 1) (B3) 

Because of the rotational invariance of H t , tj = unless m = n. So a T = Josephson 
effect can only occur for [mmm] states. This is the result of Wen and Zee. Evidently such a 
Josephson effect does not occur for the [331] state. Of course, the absence of the Josephson 
effect for the [331] state was demonstrated by different methods in sec. [Tl]. See discussion 
below eqn. ^8|. 

Several comments should be made about the above argument. The first is that one 
may readily include disorder into the argument. For example, suppose that the interwell 
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tunneling t(z) is a random function of z which fluctuates about a mean value < t(z) > 
and which is autocorrelated on some distance scale £. In this case a tunnelling event can 
change L z by an amount of order R/£ where R is the radius of the electron droplet. So the 
no-Josephson argument fails if 

(iV T -iV i )(m-n)< J R/e (B4) 

An alternative way of writing this is to note that N a = irR 2 a n a and that the difference in &_f 
of the electrons in the two edge channels is Akp = — R^/l^ so eqn. |B4| becomes 



(fc T - h) < \/i (B5) 

The interpretation of this result is simple, all the tunnelling which occurs between the two 
wells for small voltage bias (i.e. eV < 2Ai) occurs at the edges. Moreover, because of 
the requirement of momentum conservation, interedge tunneling can only occur because of 
disorder effects or because the external field is tilted. 
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FIGURES 

FIG. 1. Schematic drawing of the double quantum well. We ignore the finite well width d w . 
The right and left wells are denoted a =| and a =J_, where a is the well index (or pseudospin 
index). The external components of the field B perpendicular and parallel to the wells are denoted 
by B± and Bn . 

FIG. 2. Pair spin correlation function data from the Monte Carlo runs. The solid line through 
the circles is the analytic result for the [330] state presented in eqn. (p^). The solid line through 



the squares is the fit to the expansion given in eqn. (29) to the data. This data is for 200 particles 
per well and 2 x 10 6 Monte Carlo cycles. 

FIG. 3. The integrated tunneling conductance spectral moment Wo(-B||) associated with in- 
ter-Landau level excitations of the [330] state (solid) and the [331] state (dashed). 

FIG. 4. Schematic drawing illustrating single electron tunneling and the origin of W nc . Because 
of momentum conservation, the electron tunnels along field lines. In a tilted field this reduces the 
matrix element associated with the overlap of the single-particle Gaussian wavefunctions. This 
figure based on an argument by J. Eisenstein. 

FIG. 5. Schematic drawing of tunneling into the correlation hole of a [mmn] state with 
Fewer final states are blocked by the correlation hole than would be blocked by an uncorrelated 
electron liquid. This effect tends to increase Wq and decrease < eV >. 

FIG. 6. (a.) A plot of the radial distribution function g^(d#) vs d* = B\\d/ B±l m , the lateral 
displacement experienced by an electron during a tunneling event, (b.) A plot of correllation 
enhancement factor R vs d*. Observe that R has a minumum near the maximum of <7Tj,(d>) i.e. at 
the position of the first coordination shell. From tunneling data, one can measure R(Bn ) which, in 
general, is expected to track as it does here for the [331] state. 
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FIG. 7. Comparision of mean voltage bias results obtained for the [330] and [331] states. The 
mean voltage bias < eV > is the mean voltage of the intra-Landau level band which is observed in 
the tunneling conductance spectrum. Results are presented as a function of (d/l m )(B\\/B±) where 
d/l m = 2.4. 

FIG. 8. Mean voltage bias vs. B\\d/B±l m : (a.) Results for [330] (i.e. a pair of v = 1/3) for 
various well spacing d = 0,0.5, 1.0, 1.5, 2.4. (b.) Same as (a.) except for [550] state, (c.) Same as 
(a.) except for [331] state. 
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